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Abstract. We consider a system of N particles with a stochastic dynamics introduced by 
Brunet and Derrida [3- The particles can be interpreted as last passage times in directed 
percolation on {1, . . . , N} of mean- field type. The particles remain grouped and move like a 
traveling front, subject to discretization and driven by a random noise. As A^ increases, we 

f^ ■ obtain estimates for the speed of the front and its profile, for different laws of the driving noise. 

As shown in [7], the model with Gumbel distributed jumps has a simple structure. We establish 

pq ■ that the scaling limit is a Levy process in this case. We study other jump distributions. We 

^ I prove a result showing that the limit for large N is stable under small perturbations of the 

t3 ' Gumbel. In the opposite case of bounded jumps, a completely different behavior is found, 

where finite-size corrections are extremely small. 



m 



1. Definition of the model 



P^ . We consider the following stochastic process introduced by Brunet and Derrida [7]. It consists 

Ph ! in a fixed number A^ > 1 of particles on the real line, initially at the positions Xi{0), . . . , X]\f{0). 

^ I With {^ij{s) '■ 1 < i, j < N, s > 1} an i.i.d. family of real random variables, the positions evolve 

c3 ■ as 

X,(t + 1) = max^ {X,(t) + a,(t + 1)}. (1.1) 

^p. I The components of the A^- vector X{t) = (Xj(t),l < i < N) are not ordered. The vector 

>- ! X{t) describes the location after the t-th step of a population under reproduction, mutation 

'■^ [ and selection keeping the size constant. Given the current positions of the population, the 

^^ • next positions are a iV-sample of the maximum of the full set of previous ones evolved by an 

CN ■ independent step. It can be also viewed as long-range directed polymer in random medium 

cn ! with N sites in the transverse direction, 
O 

X,{t) = max {X^M + J2^is,is-iis); 1 < J. < iVVs = 0, . . . t - 1, j^ = z}, (1.2) 



K> I as can be checked by induction {1 < i < N). The model is long-range since the maximum in 
fll.ip ranges over all j's. For comparison with a ' 
j neighbor of i in Z in fll.ip would define the sta 
with passage time ^ on edges in two dimensions. 



H ' fll.ip ranges over all j's. For comparison with a short-range model, taking the maximum over 
■ ■ j neighbor of i in Z in fll.ip would define the standard oriented last passage percolation model 
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By the selection mechanism, the A^ particles remain grouped even when A^ — t- oo, they are 
essentially pulled by the leading ones, and the global motion is similar to a front propagation 
in reaction-diffusion equations with traveling waves. Two ingredients of major interest are: (i) 
the discretization effect of a finite A^, (ii) the presence of a random noise in the evolution. Such 
fronts are of interest, but poorly understood; see [21] for a survey from a physics perspective. 

Traveling fronts appear in mean-field models for random growth. This was discovered by 
Derrida and Spohn [13] for directed polymers in random medium on the tree, and then extended 
to other problems 



The present model was introduced by Brunet and Derrida in [7] to compute the corrections for 
large but finite system size to some continuous limit equations in front propagation. Corrections 
are due to finite size, quantization and stochastic effects. They predicted, for a large class of such 
models where the front is pulled by the farmost particles [71 [S] , that the motion and the particle 
structure have universal features, depending on just a few parameters related to the upper 
tails. Some of these predictions have been rigorously proved in specific contexts, such as the 
corrections to the speed of the Branching Random Walk (BRW) under the effect of a selection 
[1], of the solution to KPP equation with a small stochastic noise [21], or the genealogy of 
branching Brownian motions with selection [3]. For the so-called A^-BBM (branching Brownian 
motion with killing of leftmost particles to keep the population size constant and equal to N) 
the renormalized fluctuations for the position of the killing barrier converge to a Levy process 
as A^ diverges 



We mention other related references. For a continuous-time model with mutation and selection 
conserving the total mass, the empirical measure converges to a free boundary problem with a 
convolution kernel [TSj. Traveling waves are given by a Wiener- Hopf equation. For a different 
model mimicking competition between infinitely many competitors, called Indy-500, quasi- 
stationary probability measures for competing particles seen from the leading edge corresponds 
to a superposition of Poisson processes [28] . For diffusions interacting through their rank, the 
spacings are tight [27], and the self-normalized exponential converge to a Poisson-Dirichlet law 
[To] . In [1], particles jump forward at a rate depending on their relative position with respect 
to the center of mass, with a higher rate for the particle behind: convergence to a traveling 
front is proved, which is given in some cases by the Gumbel distribution. 

We now give a flavor of our results. The Gumbel law G(0, 1) has distribution function P(^ < 
x) = exp —e~^,x &M.. In [7] it is shown that an appropriate measure of the front location of 
a state X G M^ in this case is 

and that $(X(t)) is a random walk, a feature which simplifies the analysis. For an arbitrary 
distribution of ^, the speed of the front with A^ particles can be defined as the almost sure limit 

vn = lim t-^^X{t)) . 

t—^oo 

We emphasize that A^ is fixed in the previous formula, though it is sent to infinity in the next 
result. Our first result is the scaling limit as the number A^ of particles diverges. 

Theorem 1.1. Assume C,i,j{t) ~ G(0, 1). Then, for all sequences rriN —t- oo as A^ —)■ cxd, 

rriN/lnN ^ ^'^^ 
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in the Skorohod topology with S{-) a totally asymmetric Cauchy process with Levy exponent ipc 



from i3.21\j , where 

13 N = In 6jv + Nb]^^ In mjy, 
withlnhN = lniV + lnlnAr- ^ + 0(-4t7), see (E 



Fluctuations of the front location are Cauchy distributed in the large A^ limit. Keeping A^ fixed, 
the authors in [7] find that they are asymptotically Gaussian as t — )■ oo. We prove here that, 
as A^ is sent to infinity, they are stable with index 1, a fact which has been overlooked in [7|. 
When large populations are considered, this is the relevant point of view. The Cauchy limit also 
holds true in the boundary case when time is not speeded-up (mTv = 1) and N ^ oo. For most 
growth models, finding the scaling limit is notoriously difficult. In the present model, it is not 
difficult for the Gumbel distribution, but remains an open question for any other distribution. 

We next consider the case when ^ is a perturbation of the Gumbel law. Define e{x) G [— oo, 1] 
by 

e(x) = l + e^lnP(e<x). (1.3) 

Note that e = is the case of .^ ~ G(0, 1). The empirical distribution function (more precisely, 
its complement to 1) of the A^-particle system fll.ll) is the random function 

N 

UN{t,x) = N-'J2^^dt)>^ (1-4) 

This is a non-increasing step function with jumps of size 1/A^ and limits f/Ar(t, — oo) = 1, 
UN(t, +oo) = 0. It has the shape of a front wave, propagating at mean speed vn, and it 
combines two interesting aspects: randomness and discrete values. We will call it the front 
profile, and we study in the next result its relevant part, around the front location. 

Theorem 1.2. Assume that 

lim e{x) = 0, and e{x) e [-6^^,1 ~ S], (1.5) 

x—^+oo 

for all X and some 5 > 0. Then, for all initial configurations X(0) G M^, all k > 1, all 
Kn C {1, . . . , A^} with cardinality k, and all t > 2 we have 

(x,.(t)-$(X(t-l));JGA'^)^G(0,l)®^ A^^oo, (1.6) 

with $ from Ii3.20\) . and moreover, 

UN{t, ^{X{t-1)) +x) — y u{x) = 1 - e-'~^ (1.7) 

uniformly in probability as A^ —t- oo. 

As is well known, it is rare to find rigorous perturbation results from exact computations for 
such models. For example, the above mentioned, last passage oriented percolation model on the 
planar lattice, is exactly sovable for exponential passage times [IS] or geometric ones [2] on sites, 
and the fluctuations asympotically have a Tracy- Widom distribution. However, no perturbative 
result has been obtained after a decade. Even though our assumptions seem to be strong, it is 
somewhat surprising that we can prove this result. The second condition is equivalent to the 
following stochastic domination: there exist finite constants c < d {c = ln6,d = ln(l -|- S~^)) 
such that 

^ + C <sto e <sto ^ + c/, ^~G(0,1). (1.8) 

This condition is reminiscent of assumption (1.13) in [23] used to control the fluctuations of the 

front location for KPP equation in random medium. By Theorem II. 2^ as A^ — )■ oo, the front 
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remains sharp and its profile, which is defined microscopically as the empirical distribution 
function of particles, converges to the Gumbel distribution as A^ — t- oo. Hence the Gumbel 
distribution is not only stable, but it is also an attractor. 

Finally, we study the finite-size corrections to the front speed in a case when the distribution 
of ^ is quite different from the Gumbel law. 

Theorem 1.3. Let b < a and p G (0, 1), and assume that the ^i,j{t) 's are integrable and satisfy 

P(e > a) = P(e G (6, a)) = 0, P(e = a) = p, ¥{^ e{b-e, b]) > (1.9) 

for all £ > 0. Then, as N ^ oo, 

yj, = a-{a-b)il- pf'2'' + o((l - pf'2'^). 

We note that in such a case, in the leading order terms of the expansion as A^ — ;■ oo, the value 
of the speed depends only on a few features of the distribution of ^: the largest value a, its 
probability mass p and the gap a — b with second largest one. All these involve the top of the 
support of the distribution, the other details being irrelevant. Such a behavior is expected for 
pulled fronts. 

Though the mechanisms are different, we make a parallel between the model considered here, 
and the BRW with selection, in order to discuss the Brunet-Derrida correction of the front speed 
V]\f with respect to its asymptotic value. For definiteness, denote by rj the displacement variable, 
assume that i] is a.s. bounded from above by a constant a, and assume the branching is constant 
and equal to /3 > 1. The results of [1] are obtained for /3 x P(?7 = a) < 1 (Assumption A3 
together with Lemma 5 (3) in [5]), resulting in a logarithmic correction: This case corresponds 
to the Gumbel distribution for ^ in our model, e.g., to Theorems 11.11 and 11.21 . In contrast, the 
assumptions of Theorem 11.31 yield a much smaller correction (of order exponential of negative 
A^^). This other case corresponds for large A^ to the assumption /3 x P(?7 = a) > 1 for the BRW 
with selection, where the corrections are exponentially small [TT], precisely given by p^ with 
p < 1 the extinction probability of the supercritical Galton- Watson process of particles located 
at site ta at time t. In our model, the branching number is A^ and p is itself exponentially 
small, yielding the correct exponent of negative A^^, but not the factor 2^ . 

The paper is organised as follows. Section 2 contains some standard facts for the model. Section 
3 deals with the front location in the case of the Gumbel law for ^. In Section 4, we study the 
asymptotics as A^ — )■ oo of the front profile (for Gumbel law and small perturbations), and their 
relations to traveling waves and reaction-diffusion equation. In Sections 5 and 6, we expand 
the speed in the case of integer valued, bounded from above, ^'s, starting with the Bernoulli 
case. Theorems 11.11 and 11.31 are proved in Sections 13.31 and 16.31 respectively. 

2. Preliminaries for fixed A^ 

For any fixed A^, we show here the existence of large time asymptotics for the A^-particles 
system. It is convenient to shift the whole system by the position of the leading particle, 
because we show that there exists an invariant measure for the shifted process. 

The ordered process: We now consider the process X = (X(t),t G N) obtained by ordering 
the components of X(t) at each time t, i.e., the set {Xi(t),X2{t), . . . ,X]\f(t)} coincides with 
{Xi{t),X2{t), ...,XN{t)} and Xi{t) > Xiit) > ■ ■ ■ > XNif). Then, X is a Markov chain with 
state space 

A^ := {y G M^ : yi > 1/2 > . . . > Vn}- 
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Given X{t), the vector X{t) is uniformly distributed on the A^! permutations of X{t). Hence, 
it is sufficient to study X instead of X. It is easy to see that the sequence X has the same law 
as F = (Y(t),t > 0), given by as a recursive sequence 

Y{t + 1) = ordered vector^ max {Yj{t) + ^ij{t + 1)}, 1 < i < A^V (2.10) 

Note that, when X(0) is not ordered, ^(1) is not a.s. equal to ^(1) starting from Y{0) = X{0). 
In this section we study the sequence Y, which is nicer than X because of the recursion fl2.1Up : 
Denote by T^{t+i) the above mapping Y(t) ^-)■ Y(t + 1) on An, and observe first that 

F(t) = %)... T5(2)T5(i)r(0). (2.11) 

For y, X G An, write y < x ii yi < Xi for all i < N. The mapping T^(i) is monotone for the 
partial order on An, i.e., for the solutions Y, Y' of (12.101) starting from Y{0), Y'{0) we have 

F(0) < F'(0) =^ Y{t) < Y'{t), 

and moreover, with 1 = (1, 1, . . . , 1), r G M and y G M^, 

%)(y + rl) = rl + %)(y). (2.12) 

The process seen from the leading edge: For each x G M^, we consider its shift x^ by the 
maximum, 

X? = Xj — max Xj, 

and the corresponding processes X^,Y^. We call X^,Y^, the unordered process, respectively, 
the ordered process, seen from the leading edge. Note that T^(^t){y^) = T^{t){y) — (niaxj yj)l by 
ffZT^ . which yields 

a similar relation holds for x's instead of y's. Then X^, Y^ are Markov chains, with Y'^ taking 
values in A^ := {y G An : yi = 0}, and we denote by Ut the law of Y^(t). 

Proposition 2.1. There exists an unique invariant measure v for the process Y^ seen from the 
leading edge, and we have 

lim ut = u. (2.13) 

i— ^oo 

Furthermore, there exists a 5n > ^ such that 

\\ut-y\\TV<{l-SNY- (2.14) 

Similar results hold for the unordered process X°, by the remark preceeding (I2.10p . Also, we 
mention that the value of 5n is not sharp. 

Proof. Consider the random variable 

r = inf {t > 1 : ^i,i(t) = max{ei,j(t); j <N}\li<N]. 

Then, r is a stopping time for the filtration {J^t)t>o, with J^t = o-Hiji.^)', ■§ < i,i,3 > !}• It is 
geometrically distributed with parameter not smaller than 

5n = {l/Nf. (2.15) 

Denote by ©, the configuration vectors 

© = (0,0, ...,0), e = (0,-00,..., -oo). 
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They are extremal configurations in (the completion of) A^ since < |/ < © for all ?/ G A^. 
Now, by definition of r and f l2.10p . 

TgM© = T^M© = T^(r)y Vy G A^. 

Hence, for allt > r and all y G Ajv such that maxi<j<7v |/j = maxi<j<7vl^(0), 

Y{t) =T^^t)...T^^2)T^^i)y. 

We can construct a renewal structure. Define ri = r, and recursively for k > 0, t^+i = t^+toO^^ 
with 6 the time-shift. This sequence is the success time sequence in a Bernoulli process, we 
have 1 < Ti < r2 < . . . < Tfe < . . . < cxo a.s. The following observation is plain but fundamental. 

Lemma 2.1 (Renewal structure). The sequence 

(F°(s); 0<s< n), (F°(ri + s); 0<s< r^-n), (F°(r2 + s); < s < TS-T2), . . . 

is independent. Moreover, for all k > 1, (F°(rfc+s); s > 0) has the same law as (F°(l+s); s > 0) 
starting from F°(0) = ©. 

Proof, of Lemma 12. 1[ By the strong Markov property, the Markov chain Y^ starts afresh from 
the stopping times ti < T2 < . . .. This proves the first statement, and we now turn to the 
second one. Note that T^© = T^© if, for all i, {^i,j',j < A^) is a permutation of {j]ij]j < N). 
Hence, 

P(F°(1) G ■ ,ri = 1 |F°(0) = ©) = P(F°(1) G ■ |r°(0) = ©) x P(ri = 1), 

and so 

p(r°(i) G ■ |y°(o) = ©, n = 1) = p(r°(i) g ■ |f°(o) = ©) . 

From the markovian structure and by induction it follows that 

P((F°(l + s);s>0) G- |F°(0) = ©) = P((F°(l + s); s > 0) G ■ |F°(0) = ©, ri = l) 

= P((F°(l + s);s>0) G ■ |F°(0) = z, n = l) 

= P((F°(l + s);s>0) G ■ |ri = l) 

= P((F°(ri + s);s>0)G-), 

for all z G A^. I 

The lemma implies the proposition, with the law v given for a measurable F : A^ — )■ M+ by 

^ ' Tl<t<T2 



I 



^5]E(F(y°(t))w,|n = i). (2.16) 



^ ^' i>l 



Remark 2.2. (^zj T/ie proposition shows that the particles remain grouped as t increases, i.e., 
the law of the distance between extreme particles is a tight sequence under the time evolution. 
In Theorem \l.^ we will see that when the law of ^ is close to Gumhel, they remain grouped too 
as N increases. 
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(a) The location of front at time t can he described by any numerical function ^(Y(t)) or 
$(y(t — 1)) (or equivalently, any symmetric function of X(t) or X(t — 1)) which commutes to 
space translations by constant vectors, 

$(y + rl) = r + $(i/) , (2.17) 

and which is increasing for the partial order on M^. Among such, we mention also the maximum 
or the minimum value, the arithmetic mean, the median or any other order statistics, and the 
choice in li3.2(J\) below. For Proposition \2.1\ we have taken the first choice - the location of the 
rightmost particle - for simplicity. Some other choices may be more appropriate to describe the 
front, by looking in the bulk of the system rather than at the leading edge. For fixed N all such 
choices will however lead to the same value for the speed vn of the front, that we define below. 

Note that for a function $ which satisfies the commutation relation (12.1 7p we have the inequal- 
ities 

$(©) +rninmax{y,(0) +6,,(1)} < HY{1)) < $(©) + max {F,(0) + ^.(l)} • 

Now, by equation (I2.16P and by the fact that ri is stochastically smaller than a geometric 
random variable with parameter (1/A^)^ we conclude that if ^ G L^, 1^(0) G L^ then $(y(t)) G 
W, and also J \yNY'dv{y) < oo. The following corollary is a straightforward consequence of the 
above. 

Corollary 2.1 (Speed of the front). If ^ ^ L} , the following limits 

Vn = lim t~^ max{Xi(t); 1 < i < iV} = lim t"^ min{Xi(t); 1 < i < iV} 

t—^oo t—^oo 

exist a.s., and vjy is given by 

vn= dv{y) E max {yj + ^y(l)}. 

Moreover, if C, & L"^, 

t-^/2(max{Xi(t); 1 < i < iV} -VNt) 

converges in law as t ^ oo to a Gaussian r.v. with variance a^ G (0, oo). 

We call Vn the speed of the front of the X-particle system. 

Proof. The equality of the two limits in the definition of vn follows from tightness in Remark 
12. 2^ (i), and the existence is from the renewal structure. Similarly, we have 



VN= f i^HT^y) - Hy)) dv{y) 



for all $ as in Remark 12.21 (ii), where A^ is defined just before Proposition 12.11 The second 
formula is obtained by taking $(?/) = maxj<Ar?/j. The Gaussian limit is the Central Limit 
Theorem for renewal processes. I 
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3. The Gumbel distribution 

The Gumbel law G(a, A) with scaling parameter A > and location parameter a G M is defined 
by its distribution function 

P(e < x) = exp ( - e-^(^-")) , xeR. (3.18) 

This law is known to be a limit law in extreme value theory [20]. In [7], Brunet and Derrida 
considered the standard case a = 0, A = 1, to find a complete explicit solution to the model. 
In this section, we assume that the sequence ^ij is G(a, A)-distributed, for some a, A. Then, 
( = A(^-a) ~ G(0,1), while 

exp(— (^) is exponentially distributed with parameter 1, (3.19) 

and exp(— e"'') is uniform on (0, 1). Conversely, if U is uniform on (0, 1) and £ exponential of 
parameter 1, then — Alnln(l/?7) and ln£^~^ are G(0, A). 

Here, the Gumbel distribution makes the model stationary for fixed A^ and allows exact com- 
putations. 

3.1. The Front as a random walk. In this section, we fix A^ > l,a G M, A > 0. We will 
choose the function $ : M^ — )■ M, 

N 

<l>{x) = A"^ In Y^ exp Ax^ (3.20) 

to describe the front location $(X(t)) at time t. 

Theorem 3.1 ([Z])- Assume theC,i,j's are Gumbel G [a, X)- distributed. 

(i) Then, the sequence ($(X(t));t > 0) is a random walk, with increments 



^"""(H 



T = a + A-Mn V^rM (3.21) 



where the Si are i.i.d. exponential of parameter 1. 
(a) Then, 



N \ / N 



VN = a + \-^E In J2 ^r^ h (^N = A" Var In ^ Sr^ . (3.22) 



^0 ^ AO 



(Hi) The law v from proposition \2.1\ is the law of the shift V G A^ of the ordered vector V 
obtained from a N -sample from a Gumbel (7(13, \). 

Proof : Define J^t = ^(6j(s), s < t, z, j < iV), and £,,,(t) = exp{-A(e.,,(t) - a)}. By ([TI]), 

XAt + l) = max |X,(t) + a- A"Mnfi,(t + l)| 

= a + $(X(t))-A"Mn£i(t + l), (3.23) 

where 

SAt + 1) = min {S, At + i)e-^^.W|e^*(^W), t > 0. 
i<i<A' '- '■' -' 

Given J-^, each variable Si{t-\- 1) is exponentially distributed with parameter 1 by the standard 

stability property of the exponential law under independent minimum, and moreover, the whole 

vector {Si{t + l),i < N) is conditionnally independent. Therefore, this vector is independent 

of J-(, and finally, 

{Si(t), 1 < i < N,t > 1) is independent and identically distributed 
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with parameter 1, exponential law. Hence, the sequence 



N 



T(t) = a + A-Mn^fi(t)-M, t > 1, 



,j=i 



is i.i.d. with the same law as T. Now, by (13.23 p . 

$(X(t)) = ^X{t-1)) + T{t) 

t 
= $(X(0)) + 5^T(.) 

s=l 

which shows that ($(X(t));t > 0) is a random walk. Thus, we obtain both (i) and (ii). 

From (13.231) . we see that the conditional law of X{t + 1) given J-^ is the law of a TV-sample from 
a Gumbel G(a + ^{X(t)), A). Hence, i' is the law of the order statistics of a A^-sample from a 
Gumbel G(0, A), shifted by the leading edge. ■ 

We end this section with a remark. Observe that the other max-stable laws (WeibuU and 
Frechet) do not yield exact computations for our model. Hence, the special role of the Gumbel 
is not due to the stability of that law under taking the maximum of i.i.d. sample, but also to 
its behavior under shifts. 

3.2. Asymptotics for large A^. In this section we study the asymptotics as A^ — )■ oo with a 
stable limit law. When a = and A = 1, Brunet and Derrida [7] obtain the expansions 

„„ = l„7V + lnl„7V + !^ + i5^ + o(J^). (3.24) 

.^ = 5^ + .... (3.25) 

by Laplace method for an integral representation of the Laplace transform of T. We recover here 
the first terms of the expansions from the stable limit law, in the streamline of our approach. 

We start to determine the correct scaling for the jumps of the random walk. First, observe that 
S~^ belongs to the domain of normal attraction of a stable law of index 1. Indeed, the tails 
distribution is 

P(£:-^ > x) = 1 - e"^/"^ ~ x-\ x^+oo. 
Then, from e.g. Theorem 3.7.2 in [Ti] . 

where Bn = iVE(£~^;£^^^ < A^), and S is the totally asymmetric stable law of index a = 1, 
with characteristic function given for n G M by 

Ee™'^ = exp(/"(e^"^-l)^+ Ae*"^-l-znx)^ 
iJi ^ Jo ^ 

= exp < iCu — — |n||l + z— sign(u) In |m|} > 

=: exp^cH, (3.27) 
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for some real constant C defined by the above equality. By integration by parts, one can check 
that, as A^ — )■ oo, 

/"°° e~y 1 1 

hN = N —dy = N{\nN-^ + - + 0{—)), (3.28) 

Ji/N y N N^ ' 

with 7 = — /q°° e~^ \a.xdx the Euler constant. Then, 

In 67V = In iV + In In iV - -^ + C(^— ) 

IniV ^\n^N' 

We need to estimate 

^ N 

ElnV^r^-lnfe^ = Eln(l + — 5^^)) 

= mn(l + ^5)+0((J^)-), (3.29) 

for all 6 G (0, 1]: indeed, since the moments of S^^^ of order 1 — 6/2 are bounded (Lemma 5.2.2 
in [IT]), the sequence (^)^~'^[ln (l + -^S^^^^ — In (l + ^5)] is uniformly integrable, and it 
converges to 0. A simple computation shows that 



E 



ln(l + sS) = r ln(l + ey)^il + o(l)) + 0(e) ~ e ln(£-i) 

Ji y 



as £ \ 0. With e = N/h^, we recover the first 2 terms in the formula (I3.24p for v^. (If we 
could improve the error term in 03.291) to o(lnln A^/ln A^), we would get also the third term.) 
With a similar computation, we estimate as A^ — )■ 00 

Varfln(l + ^5(^))^ 
V On 



a 



N 



/ A^ N 

V„(,„(l + _5)) 

On 

00 

ln^(l + ^) 
1 



2n I y \'^y 

\nNy^ 



00 



in^ri 



y ^dy 





Co 
In AT 



In A^' 1/2 



withCo = /o°°ln^(l + y)f = 7rV3. 

3.3. Scaling limit for large A^. In this section, we let the parameters a, A of the Gumbel 
depend on A^, and get stable law and process as scaling limits for the walk: In view of the 
above, we assume in this subsection that ^ij ~ G(a, A) where a = a^ and A = \n depend on 

\ N 1 

An — T— 



bN In AT' (3.30) 

a^ = -C-\]^Hn{bN) = -C-ln2A^-(lnA^)(lnlnA^) + o(l), 
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with the constant C from fl3.27p . Correspondingly, we write 



N 



X = xW,T^(t) = a^ + A^Mn ^£:.(t)-^ • 



,i=l 

Note that, with S^^^ defined by the left-hand side of ([326]), we have by ([MI]), 

N 



j=i 
A/v V 6!v / 



\Ar V Cat 

-i^ 5o, (3.31) 

as A^ — > oo, where the stable variable Sq = S — C has characteristic function 

vr 
E exp iuSq = exp "^q^u), "^q^u) = 1 n | — m In | n | 

from the particular choice of C. In words, with an appropriate renormalization as the system 
size increases, the instantaneous jump of the front converges to a stable law. For all integer n 
and independent copies 5o,i, . . . iSo,n of Sq, we see that 

<Soi + . . . + Son , law ^ 

Inn = oq 

n 

from the characteristic function. Consider the totally asymmetric Cauchy process (iSo(r);r > 
0), i.e. the independent increment process with characteristic function 

Eexp{m(5o(r) - Sq{t'))} = exp{(r - r')^o(^)}, n G M, < r < r. 

It is a Levy process with Levy measure x"^ on M^., it is not self-similar but it is stable in a 
wide sense: for all r > 0, 

So{t) law (-, , . 

Inr = So{l) 

T 

with iSo(l) ="50. We refer to [B] for a nice account on Levy processes. 

We may speed up the time of the front propagation as well, say by a factor m^v — ^ oo when 
A^ — )■ oo, to get a continuous time description. Then, we consider another scaling, and define 
for r > 0, 

, , 1'(XW([m^r]))-$(XW(0)) 

niN 

= ^t=i ^^^^) _r\nm,, (3.32) 

rriN 

by theorem 13. 1[ Of course, this new centering can be viewed as an additional shift in the formula 
(lOnjl for ttN. By ([331]), the characteristic function xn{u) := Ee*"^^ = exp{^o(M)(l + o(l))}, 
where o(l) depends on u and tends to as A^ — )■ oo. Then, 

Kexp<.tu{ — -^ r in m^r ) > = iXNiu/mN))^ ^expj— « InmTvj 

m^v I ?TlAr 

-^ expr*o(M), 
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as A^ — !■ oo, showing convergence at a fixed time r. In fact, convergence liolds at tlie process 
level. 

Theorem 3.2. As N ^ oo, the process </?Ar(-) converges in law in the Skorohod topology to the 
totally asymmetric Cauchy process Sq{-). 

Proof. The process ^n{') itself has independent increments. The result follows from general 
results on triangular arrays of independent variables in the domain of attraction of a stable 
law, e.g. Theorems 2.1 and 3.2 in [18]. I 

Proof of Theorem li.il Apply the previous Theorem 13.21 after making the substitution C, = 

AAr(^-aAr). ■ 

4. The front profile as a traveling wave 
Recall the front profile 

N 

[/^(t,x) = iV-iJ]lx,(t)>. (4.33) 

which is a wave-like, random step function, traveling at speed v^. One can write some kind 
of Kolmogorov-Petrovsky-Piscunov equation with noise (and discrete time) governing its evo- 
lution, see (7-10) in [7] and Proposition 14. II Let F denote the distribution function of the ^'s, 
F{x) = F{^ij{t) < x). Given J^t~i, the right-hand side is, up to the factor A^~^, a binomial 
variable with parameters A^ and 

N 

P(X,(t) > x| J-i_i) = l-J]P(X,(t-l)+eij(t)<x|J-i_i) (bydnD) 

i=i 

= 1-exp-N lnF{x~y)UNit-l,dy). (4.34) 

Jr 

4.1. Gumbel case. Starting with the case of the Gumbel law F{x) = exp — e"'^'-^""-', we 
observe that fl4:M|l and fl3:20D imply 

F{Xi{t) < x|Ji_i) = exp-e^(""'^"*(^(*"^»), 

that is (13.231) . It means that X{t) — $(X(t — 1)) is independent of J-i-i, and that it is a A^- 
sample of the law G{a, A). For the process at time t centered by the front location ^{X(t — 1)), 
the product measure G{a, A)®^ is invariant. We summarize these observations: 

Proposition 4.1 (|7j). Let C,i,j{t) ~ G(a, A) be given, and X defined by U.l\} . Then, the 
random variables Gi{t) defined by G(t) = {Gi{t);i < N) and 

X{t) = G{t) + ^{X{t-l))l, t>l, (4.35) 

are i.i.d. with common law G(a, A), and G{t) is independent of X(t — l),X(t — 2),.... In 
particular, {Gi(t);i < N,t > 1) is an i.i.d. sequence with law G{a,\), independent of X{0) G 
M.^ . Moreover, 

1 ^ 
f/jv(t,x) = -^l{G,(t)>x-$(X(t-l))}, t>l,xeM. (4.36) 
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Remark 4.2. (i) The recursion lli4-36\ ) is the reaction-diffusion equation satisfied by Un- This 
equation is discrete and driven by a random noise {G{t);t > 0). 

(a) Note that the centering is given by a function of the configuration at the previous time 
t — 1. One could easily get an invariant measure with a centering depending on the current 
configuration. For instance, consider 

X{t) — inaxXj{t) = g — m.a.xgj, 
j j 

with gj i.i.d. G {a, \) -distributed, or replace the maximum value by another order statistics. 
However our centering, allowing interesting properties like the representation lj^4.!^5\ ), is the 
most natural. 

By the law of large numbers, as A^ — )■ oo, the centered front converges almost surely to a limit 
front, given by the (complement of) the distribution function of G(a, A), as we state now. 

Proposition 4.3. For allt> 1, the following holds: 

(i) Convergence of the front profile: as N ^ oo, conditionally on J^t~i, we have a.s. 

UN{t,x + <l>{X{t-l))) — > u{x) = 1 - exp(-e^^(''-")), uniformly in x G M. 
(a) Fluctuations: as N ^ oo, 

InNx {f/jv(t,x + (t-l)(ln&jv + a) + $(X(0))) - u{x)} ^ ^^{tS + tlnt + tC) 

A 

as N ^ oo, with S from \3.26\) and C from ^3.21\) . 

We wilU see in the proof that the front location alone is responsible for the fluctuations of the 
profile. It dominates a smaller Gaussian fluctuation due to the sampling. 

Proof, of Proposition 14. 3[ As mentioned above, the law of large numbers yields pointwise 
convergence in the first claim. Since U]s[{t, ■) is non-inceasing, uniformity follows from Dini's 
theorem. We now prove the fluctuation result. By fl3.2ip and (13.2 



Zn := In A^ X {$(X(t)) - $(X(0)) -tlnb^} = !^ Y^ln (l + ^S'^''\s)) 

converges in law to the sum of t independent copies of S, which has itself the law of tS+t In t+tC. 
On the other hand, we have by (14.361) . 

1 ^ 1 

UN{t + l,X + t\YibN + $(X(0))) = ^ 5^ HGi{t + 1) > X + j^^iv}. 

j=l 
By the central limit theorem for triangular arrays, for all sequences zat — )■ 0, we see that, 

1 ^ 
ArV2(_^l(Gi(t+l)>a; + ^iv}-n(a; + z^)) ^ Z ~ Ar(0, n(x)(l-n(x))) 

i=l 

as A^ ^> oo. Being of order A^~^/^, these fluctuations will vanish in front of the Cauchy ones, 
which are of order (In A^)~^. In the left hand side, we Taylor expand u{x + zn). Since G{t + 1) 
and Zjv are independent, we obtain 

InA^x {?77v(t + l,x + tln6^ + <l'(X(0))) - u{x)] -u'{x)Zn ^0 
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in probability, which proves the result. I 

Remark: A limiting react ion- diffusion equation. It is natural to look for a react ion- diffusion 
equation which has u as traveling wave (soliton). By differentiation, one checks that, for all 
V G M, u(t, x) = u{x — vt) (where u{x) = 1 — exp{— e"'^^'^'"")}) is a solution of 

ut = u^^ + A(u), (4.37) 

with reaction term 

Aiu) = A(l - u) \X In -^ + (v - A)l In -^. 

'- \ — u -' 1 — n 

Since ^(0) = A{1) = 0, the values u = and u = 1 are equilibria. For f > A, we have A{u) > 
for all u G (0, 1), hence these values are the unique equilibria u G [0, 1], with u = unstable 
and u = 1 stable. For v G [A, 3A), A is convex in the neighborhood of 0, so the equation is not 
of KPP type IISlp.2]. 

4.2. Exponential tails: front profile and traveling wave. In this section we prove Theo- 
rem [L2l We consider the case of C, with exponential upper tails, 1 — F{x) = F{^ > x) ~ e~^ as 
X — )■ +00, that can be written as 

lime(x) = 0, with e{x) = 1 + e'^lnFix). (4.38) 

(By affine transformation, we also cover the case of tails P(^ > x) ~ gA{a;-a)_-j gy definition, 
e{x) G [— oo, 1]. 

We let N -^ oo, keeping t fixed and we use $ from (I3.20p with A = 1. To show that the 
empirical distribution function (I4.33P converges, after the proper shift, to that of the Gumbel 
distribution with the same tails, we will use the stronger assumption that 

lim e{x) = 0, and e{x) G [-5~\ 1 - 5], (4.39) 

for all X with some 6 > 0. 

Proof. (TheoremOD First of all, note that \nF{x) = -(l-e(x))e-^. Let m^ = e^«(*-i)-*(^(*-i)), 
which add up to 1 by our choice of $, and let also Ei = e{x + $(X(t — 1)) — Xi{t — 1)). 
We start with the case fc = 1, K^ = {j}. From (I4.34p . 

\nF{Xj{t) - <l>(X(t - 1)) < x\J^t~i) 

N 

= J^lnF(a; + $(X(t-l))-Xi(t-l)) 

i=l 

N 

= -Y, e-"-*(^(*-^»+^'(*-i)[l - e{x + $(X(t - 1)) - Xi{t - 1))] 
1=1 

N 

1 - Y^ miSi - Y ^i^i 1 (4-40) 
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with Ii = {i : Xi(t — 1) < $(X(t — 1)) — A} and I2 the complement in {1, . . . A^}, and some 
real number A to be chosen later. By the first assumption in f l4.39p . we have 

I /J'Tijejl < sup{|£:(y)|; y>x + A}xl^O as A ^ 00, 
for fixed x. The second sum, 

will be bounded using the second assumption in fl4.39p . We can enlarge the probability space 
and couple the ,^jj(s)'s with {gij{t — l);^,j < A^), which are i.i.d. G(0, 1) independent of 
(^jj(s); i,j < N,s j^ t — 1), such that 

(7,,,(t - 1) + c < ^,j{t - 1) < (7ij(t -l) + d. 

Define for i < N, 

X,{t - 1) = m^ax{Xj(t - 2) + gij{t - 1)}. 

By the previous double inequality, 

X,{t - 1) + c < X,{t - 1) < X,{t - 1) + d, 
and, since $ is non- decreasing and such that $(y + rl) = $(y) +r, we have also 
^{X{t - 1)) - <l>(X(t - 2)) > <|)(X(t - 1)) - ^X{t - 2)) + c. 



On the other hand, in analogy to the proof of Proposition 14.11 for the Gumbel case we know 
that {Xi{t - 1) - $(X(t - 2)); 1 < z < A^) is a A^-sample of the law G(0, 1). So, 

$(X(t - 1)) - $(X(t - 2)) = ln(M+ln('l + ^5(^)') 

= lniV + lnlnA^ + o(l) 

in probability from (13. 26 p . and 

max{Xj(t -l)]i < N} - $(X(t - 2)) - In A^ converges in law 

by the limit law for the maximum of i.i.d.r.v.'s with exponential tails [201 Sect. 1.6]. Combining 
these, we obtain, as A^ — )■ 00, 

$(A:(t-l))-max{Xi(t-l);i< AT} > $(X(t - 1)) - rf - max{Xi(t - 1); i < A^} + c 

— !■ +00 in probability , 

which implies that the set I2 becomes empty for fixed A and increasing A^. This shows that 
Y^i£h e^'(*"^)"*(^(*"^)) -)■ in probability (i.e., under P(-| J't„2)) uniformly on X{t-2). Letting 
N ^ 00 and A -^ +00 in (I4.40p . we have 

¥{Xj{t) - $(X(t - 1)) < x\Tt^2) -^ exp -e"^ 

as A^ — )■ 00 uniformly on X(t — 2), which implies the first claim for k = 1. For k > 2, recall that, 
conditionally on J-t_i, the variables {Xi{t);i < N) are independent. The previous arguments 
apply, yielding ([L6D- 

Statement (ll.7p for fixed x follows from this and the fact that Xi{t) are independent condi- 
tionally on J^t-i- Convergence uniform for x in compacts follows from pointwise convergence of 
monotone functions to a continuous limit (Dini's theorem). Uniform convergence on M comes 
from the additional property that these functions are bounded by 1. ■ 
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Remark 4.4. (i) From the stochastic comparison i\1.8\} of ^ and the Gumhel, we obviously have 
Vn = In^AT + 0{1). We believe, but could not prove, that the error term is in fact o(l). 
(a) We believe, but could not prove, that the conclusions of Theorem \1.2\ hold under the only 
assumption that the function e from ^.38 ) tends to at +oo. 



5. Front speed for the Bernoulli distribution 
In this section we consider the case of a Bernoulh distribution for the .^'s, 

niiM = 1) = p, m^At) = 0) = g = 1 - p, 

with p G (0,1). For all starting configuration, from the coupling argument in the proof of 
proposition 12. 11 we see that all A^ particles meet at a same location at a geometric time, and, at 
all later times, they share the location of the leading one, or they lye at a unit distance behind 
the leading one. We set $(x) = maxjxj; j < A^}, and we reduce the process X^ to a simpler 
one given by considering 

Z{t) = tt{j : 1 < J < N,Xjit) = 1 + max{X,(t - 1);^ < N}}. (5.41) 

Z{t) is equal to the number of leaders if the front has moved one step forward at time t, and 
to if the front stays at the same location. Here, we define the front location as the rightmost 
occupied site $(X(t)) = max{Xj(t); j < N}. Then, it is easy to see that Z is a Markov chain 
on {0, 1, . . . , N} with transitions given by the binomial distributions 

P(z(..i)^,zM^™).{«[;^-;:^:|(:j;:^j; (5,42, 

Note that the chain has the same law on the finite set {1,2,...} when starting from or from 
N. Clearly, w^v — )■ 1 as iV — )■ oo. We prove that the convergence is extremely fast. 

Theorem 5.1. In the Bernoulli case, we have 

t;Ar = l-g^'2^ + o(g^'2^) (5.43) 

as N ^ oo. 

Proof. The visits at of the chain Z are the times when the front fails to move one step. Thus, 

t 
$(X(t)) = $(X(0)) + 5^1z(.)^o, 

s=l 

which implies by dividing by t and letting t — )■ oo, that 

Vn = i^n{Z 7^ 0) = 1 - ijn{Z = 0), 

where i^at denotes the invariant (ergodic) distribution of the chain Z. Let Ej^j, P^ refer to the 
chain starting at A^, and T^ = inf{t > 1 : Z{t) = k} the time of first visit a.t k {0 < k < N). 
By Kac's lemma, we can express the invariant distribution, and get: 

^^ = 1 - (^oTo)-^ = 1 - {EmTo)-\ (5.44) 

Let ctq = 0, and cri, o"2 . . . the successive passage times of Z at A^, and J\f = X]i>o '^(^i<To the 
number of visits at A^ before hitting 0. Note that A/" has a geometric law with expectation 
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EnM = Pn{Tq<TnY\ Then, 
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EnTq 



E 



j>l 

^ En [{(Ti - cri„i)l^^<rj + En{To - cr^) 
«>i 

y^ -^Af [la,_i<To-EAr(ail^i<To)] + -EAr(To|To < Tat) (Markov property) 



i>l 



EN[^f] X E^(ail,,<To) + ^7v(To|To < Tj 
1 - P^(To < Tn 



Pn{To < Tn] 
We will prove a Lemma. 
Lemma 5.1. We have 

as N tends to oo. Moreover, 



En{Tn\Tn < Tq) + _E'Ar(ro|ro < T, 



N) 



(5.45) 



N^nN 



Pn{To < Tn) ~ g^^ 2 



lim E^(To|ro<Tjv) = 2, 

Af— >oo 



lim En{Tn\Tn < To) = I. 



(5.46) 



(5.47) 



(5.48) 



The lemma has a flavor of Markov chains with rare transitions considered in [9], except for 
the state space which is getting here larger and larger in the asymptotics. With the lemma at 
hand, we conclude that 



EnTq 



1 



Pn{To < Tn) 
1 



as A^ tends to oo. From (I5.44p . this implies the statement of the theorem. 



Proof, of lemma 15.11 We start to prove the key relation (I5.46p . We decompose the event 
{To < T/v} according to the number i of steps to reach from state A^, 



P^(To < Tn) = 5^P^(To = i < Tn). 



(5.49) 



e>i 



We directly compute the contribution oi i = 1: By 05.4ip . we have 

P;v(To = l<TAr)=g^', 



(5.50) 
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which is neglegible in front the right-hand side of fl5.46p . We compute now the contribution of 
strategies in two steps: 

Af-l 

fc=i 



.A^ 



^-1 / ,r 

2 v^ fN 



k=l 



r ElV)(i-? 



k 



N\k 



g^'2^. 



(5.51) 



For £ > 2 we write, with the convention that /cq = N , 

N-l 
k-i,...ki=l 



PM{To = i+l<T, 



N 



1 < TM,Z(i) = ki,i = 1, 



I..— 1 ,— 1 V / 



ki,...ki=l i=l 

N-l e 

s E [n 

A;i,...fc£=l i=l 
N-l e 



Q 



ki.i{N-ki) 



■rkiN 



Q 



Ar2 



E n 

fci ,...kp=l i=l 



N 



Q 



ki{N-ki^i) 



(since ko = N) 



Ar2 

q ae, 



(5.52) 



which serves as definition of a^ = a£{N). For e G (0, 1), define also bt = bi{e, N) by 



E n 



l<ki,...ke<N-l,ki>(l-£)N i=l 



N 



Q 



ki{N-ki^i) 



Then, by summing over ki 



£-1 



E [n 

l<fcl,...fcf_l<Af-l 4 = 1 

e-1 



N 



Q 



fci(Af-fci_i) 



ki{N-ki^i) 



with 



^ E [nifci" 

l<fci,...A:£_i<Af-l 4=1 

= E + E 

ke-i<{l-£)N ki_i>{l-e)N 

< 7Ar a£_i + (1 + g)^6^_i. 



def 



[(l + g^-^^-i)^-l] 



(5.53) 



Ne\N 



lN = lN{e) = {l + q''T-^-Nq 



-.Ne 
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as N -^ oo. We now bound bi in a similarly manner. First we note that, for any rj such that 
T] > —e\n{e) — (1 — e) ln(l — £:) > 0, we have 

y I 7 ) < exp Nr] for large A^. 

ke>{l-e)N ^ ^ 

Note also that we can make rj small by choosing e small. Then, 



n 



Q 



ki{N-ki^i) 



^Nr,^{l-e)N{N-ki_i) 



h < E 

l<fci,...fcf_i<Af-l li=l 

= E + E 

ke-i<{l-£)N ke_T_>{l-£)N 

< g^(i-^)A^^e^''a,_i + g(i-^)^e^''6,_i. (5.54) 

For vectors u, v, we write u < v ii the inequality holds coordinatewise. In view of fl5.53p and 
(15.541) . we finally have 

where the matrix M is positive and given by 

7^ (1 + q^ 

e(l-e)N'^ ^Nrj {l-e)N^Nr, 

We easily check that, for e and rj small, M has positive, real eigenvalues, and the largest one 
A+ = X+{N,e,r]) is such that A+ ~ 7Ar as A^ — )■ oo. By (15.551) . 

ae < X^+'^{a2 + b2), 
and since A+ = X+{N,e,r]) < 1 for large A^, 



M 






0-2 + &2 



e>2 



A, 



(5.56) 



Now, we estimate 02 = 02 and 62? both of which depend on A^: 



a2 



E 



l<A;i,fc2<Af-l 



:)(:)''"^-'- 



< 



E (^) [d + «"-'=•)" -1] 



= E + E 

ki<(l-e)N ki>{l-e)N 

< 7iv2^ + (l + g)^e^^ 



and 



6, < ,<'-'>" 



E 



l<fci,fc2<Af-l,fc2>(l-e)Af 



A^\ fN 



ki I \k 



< 2^q(^-')^e^\ 



From (I5.56p . we see that 



5^a,(Ar) = o(2^), 



£>2 



and, together with (I332D . I^M), flH3U]) . (I53TD . it implies ([230]) • 
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The limit f l5.47p directly follows from the above estimates. 
Finally, we turn to the proof of (15.481) . Note that 



N 



l-il-pf) +Ejv(T^1t^.>i) . (5.57) 

We only need to show that the last term is exponentially small. For that, we use Markov 
property at time 1, 

E^(T^1t^>i) < [l- (l-(l-p)^)'^](l + maxE^(Tjv)), 

where the first factor is exponentially small. To show that the second factor is bounded, one 
can repeat the proof of part a) with p = 1 of the forthcoming Lemma 16.51 ■ 

6. The case of variables taking a countable number of values 

In this section we consider the case of a random variable C, taking the values N^ := {/ G Z : / < 
k}, with fc G Z, so that 

m^At) = I) = Pi. (6.58) 

for / G Nk, with pi > 0,Pk ^ (0, 1) and X]«eN Pi ~ ^- ^^ ^^ ^^^ Bernoulli case we can reduce 
the process X^ to a simpler one given by Z{t) := {Zi{t) : / G Nk), where 

Zi{t) = tt{j : 1 < J < N,Xj{t) = max{X,(t - 1); 1 < z < iV} + /}, (6.59) 

for / G Nfc. Note that Zk{t) is equal to the number of leaders if the front has moved k steps 
forward at time t, and to if the front moved less than k steps. Z is a Markov chain on the set 



|mG{0,...,iVf- J] 



fii. := < m G 10 iVV^'= : 2^ m,- = N 

JSNfe 

where m, are the coordinates of m. We now proceed to compute the transition probabilities of 
the Markov chain Z. Assume that at some time t we have Zt = m = {rrii : i G N^). For each 
i G Nfc, this corresponds to rrii particles at position i. Let us now move each particle to the 
right adding independently a random variable with law ,^o.o- We will assume that nik > 1. The 
probability that at time t + 1 there is some particle at position k is 

Sk{m) := 1 - I ^ Pi 

\l=—oo 

Similarly, the probability that the rightmost particle at time t + 1 is at position k — 1 is 

/ fe-l \ '"fe / fc-l \ "fc-i / fc-2 \ *"* 

Sfe_i(m) ■■= i^ Pij - {J2 Pn [J2p^ 

\l=—oo / \l=—oo / \l=—oo 

In general, for r ^N^-, the probability that at time t + 1 the rightmost particle is at position r 
is 



Srirn) 
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Define now on Qk the shift 9m by {9m)i = rrii^i for i E Nk- For r G N^, let Sr (m) := Sr{Om) 
and in general for j > 1 let 



sl,^'{m) := Sr{0^m). 



Si) I 



Define s(rra) := {sr{m) : r G N^), for j > 1, s^^\m) := {sy{m) : r G N^). Dropping the 
dependence on m of s^, Sr , s and s'^-'^ we can now write the transition probabilities of the 
process Z(t) as 



P(Z(t + l) = n\Z{t) =m) 



A^(iV;s(2))(n), mfc_2 > l,mfc = mfc_i = 0, 
A^ (iV; s^^)) (ra), mfc_j > 1, m^ = mfc_i = ■ ■ 



(6.60) 



rrik-j+i = 0, 



where for Ui with Yli'^i ~ -'-' A^(^;w) denotes the multinomial distribution (with infinitely 
many classes). Let us introduce the following notation 

k—i 
j=-oo 

for integer i > 1. 

Assumption (R). We say that a random variable ^ distributed according to (16.581) satisfies 
assumption (R) if 



Pk X pk-i > 



and 



IE(|6,o|) < oo. 
We can now state the main result of this section. 



Theorem 6.1. Let ^ be distributed according to ^6.58\) and suppose that it satisfies assumption 
(R). Then, we have that 

Vr, = k-qf2^ + o{qf2'^), (6.61) 

as N ^ oo, where qk '■= 1 — Pk- 



6.1. Proof of Theorem 16. IL To prove Theorem 16.11 we will follow a strategy similar to the 

one used in the Bernoulli case. Let us first define for each m = {irti : i G N^) G Vtk the function 

= 0(m) := sup{i G Nfc : mi > 0}. (6.62) 

As in the Bernoulli case, we denote by v^ the invariant (ergodic) distribution of the chain Z. 

Lemma 6.1. Let C, be distributed according to 116. 5 8\} . Then, we have that 



Vn = k - Un{(J) < k - I) -"Y z/7v((/) <k-j). 

i=2 



(6.63) 
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Proof. Let $(x) = max{xj; i < N}, and note that for every positive integer time t 

t 
$(X(t)) = <|.(X(O)) + 5^0(Z(^)). 

Hence 

oo 
WTV = ^ iz/Ar(0 = i) = /C - ^ Z/Ar(0 < /C - j). (6.64) 



We will now show that the first two terms of the expression for the velocity f l6.63p given in 
Lemma 16. Ij dominate the others. 

Lemma 6.2. Let C, be distributed according to 116. 5 8\} . Then, for each i >2 we have that 

z/Ar(0< A;-i) < ( — j z>7v(0< A;-l). 
Proof. Let us fix m G ^2^. Define k, := sup{i G N^ : mj > 0}. Let us first note that 

Pm{<f>iZ{l)) <k-l) = r^-^^, 
while 

P^(0(Z(1)) <k-2)= r;"«-^rr''^ < M j Pm{<f>{Z{l)) <k-l). 
Hence 

PnMZ{t)) <k-2) = Yl Prn{Z{t - 1) = TTl' , 0(Z(t)) < k ~ 2) 

= Yl Prn{Z{t - 1) = m')P^\ct>{Z{l)) <k-2) 

< [^-^] PU<PiZit))<k-l), 

where in the last inequality we used the fact that by definition m'^, > 1. A similar reasoning 
shows that in general, for i > 2, 

Pm{(t^{Z{t)) <k-l)<r-^] PM{Z{t)) <k-l). 



Taking the limit when t — )■ oo and using Proposition 12. H we conclude the proof. 
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Lemma 6.3. Let ^ be distributed according to /i6.58\) and suppose that assumption (R) is 
satisfied. Then 

Proof. Note that by summation by parts, assumption (R) implies that 

oo 

y^ji < oo. 

j=2 

Therefore, 

j^ \rij n \rij ^ y\ri^ 

■ 

Theorem 16.11 now follows from Lemmas I6.1[ I6.2[ 16.31 and the next proposition, whose proof we 
defer to subsection 16.21 



Proposition 6.2. We have that 



Vn{4> <k-l) 
lim r^ = 1. 



6.2. Proof of Proposition [6T2l Let us introduce for each m & Qk the stopping time 

T^ := inf{t > 1 : Z{t) = m}. 

Define now Q^ := {m G Qk '■ i^k = 0}. Furthermore, we denote in this section ©:=(..., 0, A^) G 
Qk- We now note that by Kac's formula 

i^N{Zk = 0) = ^ un{Z = n) = ^ -— 



neno neno ^^ ^> 



Hence we have to show that 



lim ^"^;g t:^^ = I. (6.65) 

We will prove f l6.65p through the following three lemmas. 

Lemma 6.4. Assume that ^ is distributed according to i\6.58\} . Then, for every n & VlP^ we have 
that 

E^{T^) = E^iT^.T^ < T„)— — i— — + E^{T^\T^ < T^) + Um^u), 

where 1 — e"*"^ < inf„^^...^„j._j^o I^a^I < sup„^^ „^_^ q I^^a^I ^ 2 + e~^^ for some constant C > 0. 

Lemma 6.5. Assume that ^ is distributed according to Ii6.58\) . Then, there is a constant C > 
such that the following are satisfied. 
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a) For p = 1 and p = 2, and for every N >2 we have that 

sup E„ra< 2^(1 + e-^^). (6.66) 

b) For every N > 2 we have that 



sup \E^{T^,T^ < T„) - 1| < e-^^. 
To state the third lemma, we need to define the first hitting time of the set Q^. We let 



Ta := inf T^. 

meno 

Lemma 6.6. Assume that ^ is distributed according to Ii6.58\) . Then, there is a constant C > 
such that 

J2 PeiTn < Te) = P^iTA < Te)(l + 0(6"^^)) . 

Let us now see how Lemmas 16 ■4[ 16.51 and 16.61 imply Proposition 16.21 We will see that in fact, 
Proposition 16.21 will follow as a corollary of the corresponding result for the Bernoulli case with 
q = Qk- Note that Lemma [6.41 and part (6) of Lemma [6.51 imply that 

1 - e"^^ 
Pe(T„ < Te) > p .^. , n G ^l 

Hence, summing up over n G fi^, by Lemma [6.61 we get that, for some C > 0, 

P^Ta < Te) > (1 - e-^'^) J^ -1-. (6.67) 

Now, note that Te(^A < Te) is equal to the probability to hit before A^, starting from A^, 
for the chain Z defined through random variables with Bernoulli increments as in Section 5. 
Hence, by 05.461) of Lemma [5. II we conclude that for A^ large enough 

(l + e-^%r2^>5:-i-. (6.68) 

On the other hand, applying the Cauchy-Schwarz inequality to the expectation EQ^{-\Tn < T^) 
in Lemma 16.41 and using Lemma 16. 5[ we obtain for each n ^ Q^ that 

£<- + -= + ., 



where Oi := 1 + e"^^, 02 := 2(1 + e"^^) and Og := Un, E := K(T„), P := P(s{Tn < Te) and 
we have used 06.66P of part (a) of Lemma [6.51 with p = 2. It follows that 



1 > V '^2 - 4ai(a3 - E) -a2 



/P - 2ai 

Hence, 



ai^>E-^^^/al-Aa^ias-E). 
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Now, al - 4ai(a3 - E) < 8(1 + E) for large A^, so that 



"'?^^^-^7ff (l^^ 



Now, by inequality f l6.68p we conclude that for A^ large enough -^ < q^ 2^+^. Therefore, 
Summing up over n G fi^, by Lemma [6.61 we get that 



P^{Ta < Te) < (1 + e-^'^) J^ -^^ (6.69) 



neQl 



for some C" > 0. Finally, 05.46P of Lemma 15.1^ together with inequalities (I6.67P and (I6.69p . 
imply inequality (I6.65p . which finishes the proof of Proposition 16.21 I 



6.2.1. Proof of Lemma \6. 5\ Part (a). We will first prove that there exists a constant C > 
such that 

sup Pm{T^ > 2) < e"^^. (6.70) 

The strategy to prove this bound will be to show that with a high probability, after one step 
there are at least ^^ leaders. This gives a high probability of then having A^ leaders in the 
second step. Consider now the set Lk^N '■= \j^ ^ ^k ■ f^k > [^^] }• We have 

Pm{T^<2)>p(x>^-^\ inf P„(Te = l), (6.71) 

where X is a random variable with a binomial distribution of parameters pk and A^. Now, by a 
large deviation estimate, the first factor of (I6.7ip is bounded from below by 1 — e"*"^. On the 
other hand, we have for m G Lk^N, 



P™(Te = l)>(l-(l-Pfc)^^''/') >1- 



e 



for some constant C > 0. This estimate combined with (16.711) proves inequality (I6.70p . Now, 
by the Markov property, we get that, for all m G ^2^, 

EmiTco) = -E'm(T®lTgj<2) + 2_^ -E'm(^elTe>2,Z(2)=n) 

< 2P„,(re < 2) + 5^ ^r„(lTe>2,z(2)=n[2 + P„(re)]) 



< 2P„(Te < 2) + 2 + sup EniT^) P^iT^ > 2), 

where the supremum is finite, in fact smaller than 5^^ with 6j\[ from (12.150 . Bounding the first 
term of the right-hand side of the above inequality by 2, taking the supremum over m & Qk 
and applying the bound (I6.70p . we obtain (16.661) of (a) of Lemma [6.51 with p = 1. The proof of 
(I6.66P when p = 2 is analogous via an application of the case p = 1. 
Part (b). Note that for every state ?7i G f2^ we have that 
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EeiT^^T^ < Ta) < -E'e(T®,T® < T^) < -E©(T®). 
Hence, it is enough to prove that 

|i?®(Tg5)-l|<e-^^, (6.72) 

and that 

\E^{T^,T^<TA)-l\<e-^''. (6.73) 

To prove (I6.72p note that 

E^{T^) = (1 - {l-p,ff + E^{T^,T^ > 1). (6.74) 

But by the Markov property, 

Ee(Te,Te>l)< (l-{l-{l-pkff) (l+ supE„(Tg,)). 
Note that 






AT' 



i_iV(i-p,; 

Using part (a) just proven of this Lemma, we conclude that 

i^e(Te,Te>l)<e-^^. (6.75) 

Substituting this back into (I6.74p we obtain inequahty (I6.72p . To prove inequahty (I6.73p . as 
before, observe that 

E^{T^,T^ < Ta) = (1 - {l-Pkff + E^{T^,Ta > T^ > 1). (6.76) 

Noting that E^{Tq,,Ta > T® > 1) < Eq){T^,T^ > 1), we can use the estimate (16.751) to obtain 



6.2.2. Proof of Lemma 6^. We will use the following relation, which proof is similar to that of 
(15.451) and will be not be repeated here: for every n & Q^, 

Ee(T.) = i?e(Te|Te < T„)§^j^^^|^ + i?e(r„|T„ < T^). (6.77) 

Let us now derive Lemma [6. 4[ Let n G ^2^ and m G fi^. We first make the decomposition 

EUTn) = {T)i + (T)2, (6.78) 

where 

(T)i := £'„(T„lTe<T„) and 
(T)2 := -Em(T„lrg3>r„)- 



We also denote by (r)2 the supremum of (T)2 over all possible n E fl^ and m eVL^. Now, 

(T)2 = (T)2i + (T)22, (6.79) 

where 
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(r)2i := -E'm(T„lre>r„lzfc(i)>c7v) and 

(r)22 := -E'm(^nlTffi>T„lZfc(l)<CA')- 

Now note that for any constant C < p^, by the Markov property and a standard large deviation 
estimate we have that 

(T)22 = Pm{Tn = l)+ 2j -^m(^nlT®>T„>2lZfc(l)=2i ) 

zi<CN 

< PUZ{1) = n) + (l + (T);) PmiZki^) < CN, Z{1) ^ n) 

< (l+jT)t)PUZk{l)<CN) 

< (l+Ir)l)e-^^, (6.80) 

for some constant c > depending on C,pk- On the other hand, by definition of the event 
{T® > T„,}, we have the first equahty below: 

(T)2l = -Em(T„lT(B>T„lZfe(l)>CA'lZfe(2)<Af-l) 



< (l-(l-(l-p,f^)^)(2+(T)2) 

< C'N{l-p,f''{2 + jTY,), (6.81) 

for some C > 0. We can now conclude from (16. 79 p . (I6.80p and (I6.8ip . that there is a constant 
C > such that 



(r)2 < Ce"^^. 

Let us now take m = n & Q^ and examine the first term of the decomposition (16.781) . Note 
that by the strong Markov property, 

(T)i = Ee(Tn) + E^{T^1t^^tJ. (6.82) 

Now, by part (a) Lemma l63] with p = 1, we see that the second term in the above decomposition 
is bounded above as follows, 

E„,(Te)<2(l + e-^^). (6.83) 

Collecting our estimates, we get 

EfiTn = En[Tn] T^ < T„,) + En[Tn] Tn < T^) 

= En\T(^] Tq < Tn) + En[Tn — Tq] Tq < Tn) + En[Tn', Tn < T^) 

= E„(Te; Te < Tn) + Pn{Te < Tn) x E^{Tn) + EniTn, Tn < Te). 



Here we bound the first term with (I6.83p . the last one by (T)2, and we can use (16.771) to obtain 
the desired conclusion. ■ 

6.2.3. Proof of Lemma 1 6. (A First note that 

Y, P^iTn < Te) > P^{Ta < Te), 
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and it suffices to prove an inequality in the converse direction. It is natural to introduce the 
number Ma of visits of the chain to the set fi^ before reaching the © state, 

since we have, for all m & Q^, the relations 

nenl 
Then, by the strong Markov property, 

E^i^fA) = E^{MaIna>i) 

< I 1 + sup EniAfA) I Pe(A/'A > 1) • (6.85) 

V "sno J 

In view of fl6.84|) . where the ffist term is smaller than the last one, it suffices to show that 



sup EniAfA) = 0{e 

nGQO 



~C'N\ 



in order to conclude the proof of the Lemma. In this purpose, use the strong Markov property 
to write 

EniAfA) = Er, {AfAlT^=l) + En {ArAlT^>2) 

= 0+ 2_^ En (lTA<Tg3,Z(TA)=m(l + E„,J\fA)) 

< f 1 + sup EUAfA)] PniTA < Te) . (6.86) 

Observe also that, for all n G f2fc, 

Pn{TA<T^) < P„(Ta = 1) + P„(Te > 2) 

< (1 - pfc)^+ sup P„(Te> 2) 

neflk 

< 2e-^^ (6.87) 

by (16. 701) . Now, the desired result follows from (16.861) and (16.871) . provided that the supremum 
in the former estimate is finite. To show this, note that sup^ Pm(T^ > 2) < (1 — Pk)^ , 
which implies that T^ is stochastically dominated by a geometric variable with this parameter. 
Therefore, 

SUpEm{AfA) < sup Em{T(s) < (1 - Pk)~^ , 
m m 

ending the proof. ■ 
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6.3. Proof of Theorem 11.31 Changing ^ into (^ — a)/{a — b), we can restrict to the case 
a = 0,b = —1. Then, for fixed e > 0, we define i.i.d. sequences ^ij(t) and ^ij(t) by 

Clearly, these variables are integrable since ^ is. Since ^ijit) < C,ij(t) < ^i,j{t), the correspond- 
ing speeds are such that 

vn <vn < vn ■ 
From Theorem 16.11 both vn and (1 + e)~^VN are —(1 — p)^ 2^ + o((l — p)^ 2^) as A^ — > oo, 
which, in addition to the previous inequalities, yields 

-{1+e) < liminf 17^(1 -p)-^'2-^ < limsupt;^(l - p)-^'2^^ < -1. 
Letting e \ 0, we obtain the desired claim. ■ 



References 

[I] M. Balazs, M. Z. Racz, B. Toth: Modeling flocks and prices: jumping particles with an attractive interaction. 
|http : //arxiv . org/abs/1107 . 3289| to appear in Ann. Inst. H. Poincare - Probabilites et Statistiques 

[2] J. Baik, P. Deift, Percy, K. Johansson. On the distribution of the length of the longest increasing subsequence 

of random permutations. J. Amer. Math. Soc. 12 (1999), 1119-1178. 
[3] J. Berestycki, N. Berestycki, J. Schweinsberg: Survival of near-critical branching Brownian motion. J. Stat. 

Phys. 143 (2011), 833-854 
[4] J. Berard, J.-B. Gouere: Brunet-Derrida behavior of branching-selection particle systems on the line. Comm. 

Math. Phys. 298 (2010), 323-342. 
[5] J. Berard, J.-B. Gouere: Survival probability of the branching random walk killed below a linear boundary. 

Electron. J. Probab. 16 (2011), 396-418. 

[6] J. Bertoin: Levy processes. Cambridge Tracts in Math. 121. Cambridge Univ. Press, Cambridge, 1996 
[7] E. Brunet and B. Derrida: Exactly soluble noisy traveling-wave equation appearing in the problem of 

directed polymers in a random medium. Phys. Rev. E 70 (2004), 016106 
[8] E. Brunet, B. Derrida, A. Mueller, S. Munier: Noisy traveling waves: effect of selection on genealogies. 

Europhys. Lett. 76 (2006), no. 1, 1-7. 
[9] O. Catoni, R. Cerf: The exit path of a Markov chain with rare transitions. ESAIM Probab. Statist. 1 

(1995/97), 95-144 
[10] S. Chatterjee, S. Pal: A phase transition behavior for Brownian motions interacting through their ranks. 

Probab. Theory Related Fields 147 (2010) 123-159 

[II] O. Couronne, L. Gerin: Survival Time of a Censored Supercritical Galton- Watson Process. 
|http://fr. arxiv.org/abs/llll . 1078 , to appear in Ann. Inst. H. Poincare - Probabilites et Statistiques 



[12] A. Dembo, O. Zeitouni: Large deviations techniques and applications. Applications of Mathematics 38. 

Springer- Verlag, New York, 1998 
[13] B. Derrida, H. Spohn : Polymers on disordered trees, spin glasses, and traveling waves. J. Statist. Phys., 

51 (1988), 817-840. 
[14] R. Durrett: Probability: theory and examples. Fourth edition. Cambridge University Press, Cambridge, 

2010. 
[15] R. Durrett, D. Remenik: Brunet-Derrida particle systems, free boundary problems and Wiener-Hopf equa- 
tions. To appear in Annals of Probability, 2011 
[16] B. Gilding, R. Kersner: Travelling waves in nonlinear diffusion-convection reaction. Progress in Nonlinear 

Differential Equations and their Applications, 60. Birkhauscr Verlag, Basel, 2004 
[17] I. Ibragimov, Y. Linnik: Independent and stationary sequences of random variables. Wolters-Noordhoff 

Publishing, Groningen, 1971. 
[18] J. Jacod. Thcorcmcs liniite pour les processus. Saint-Flour notes 1983, 298-409, Lecture Notes in Math., 

1117, Springer. 
[19] K. Johansson. Shape fluctuations and random matrices. Comm. Math. Phys. 209 (2000), 437-476 



LAST PASSAGE PERCOLATION AND TRAVELING FRONTS 30 

[20] M. Leadbetter, G. Lindgren, H. Rootzen: Extremes and related properties of random sequences and pro- 
cesses. Springer Series in Statistics. Springer- Verlag, New York-Berlin, 1983 
[21] P. Maillard: Branching Brownian motion with selection of the N right-most particles: An approximate 

model, http://arxiv.org/abs/1112.0266, to appear in Ann. Inst. H. Poincare - Probabilites et Statistiques 
[22] S.N. Majumdar, P. Krapivsky: Extreme value statistics and traveling fronts: Application to computer 

science Phys. Rev. E 65 (2002), 036127 
[23] C. Monthus, T. Garcl: Anderson transition on the Cayley tree as a traveling wave critical point for various 

probability distributions. J. Phys. A 42 (2009), 075002 
[24] C. Mueller, L. Mytnik, J. Quastel: Effect of noise on front propagation in reaction-diffusion equations of 

KPP type. Invent. Math. 184 (2011), 405-453 
[25] J. Nolen: A central limit theorem for pulled fronts in a random medium. Netw. Heterog. Media 6 (2011), 

167-194 
[26] D. Panja. Effects of Fluctuations on Propagating Fronts, Physics Reports 393 (2004), 87-174. 
[27] S. Pal, J. Pitman. One-dimensional Brownian particle systems with rank dependent drifts. Ann. Appl. 

Probability 18 (2008), 2179-2207 
[28] A. Ruzmaikina and M. Aizenman: Characterization of invariant measures at the leading edge for competing 

particle systems Ann. Probab. 33 (2005), 82-113 

(Francis Comets) Universite Paris Diderot - Paris 7, Mathematiques, case 7012, F-75 205 Paris 

Cedex 13, France 

E-mail address: comets@math.univ-paris-diderot.fr 

(Jeremy Quastel) Departments of Mathematics and Statistics, University of Toronto, 40 St. 
George Street, Toronto, Ontario M5S 1L2, Canada 
E-mail address: quasteKSmath. toronto . edu 

(Alejandro F. Ramirez) Facultad de Matematicas, Pontificia Universidad Catolica de Chile, 
Vicuna Mackenna 4860, Macul, Santiago, Chile 
E-mail address: ararniirezOmat . puc . cl 



